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Discussion of "Statistical Modeling of 
Spatial Extremes" by A. C. Davison, 
S. A. Padoan and M. Ribatet 



Benjamin Shaby and Brian J. Reich 



The review paper on spatial extremes by Davison, 
Padoan and Ribatet is a most welcome contribu- 
tion. The authors cover quite a lot of ground, mak- 
ing connections between different approaches while 
highlighting important differences. In particular, we 
applaud their careful attention to model checking, 
which can be difficult in general but particularly so 
for spatial extreme value models. 

1. PREDICTION AND MODEL VALIDATION 

With its extensive set of diagnostics for evaluat- 
ing model fit, this paper provides a nice template for 
practitioners to follow. However, perhaps the most 
important feature of spatial models is their ability to 
predict at unobserved sites. The account presented 
here does not address the prediction problem, which 
is both a critical task in its own right and a tool for 
comparing models. More traditional spatial analy- 
ses typically include various performance metrics to 
evaluate prediction at a withheld test set of obser- 
vation locations. While spatial prediction is difficult 
for the max-stable process models described in the 
paper, computational tools to accomplish this task 
do exist (Wang and Stoev (2011)). Spatial predic- 
tion for copula models is considerably more straight- 
forward. 

However, we note that even with predictions at 
hold-out locations in hand, evaluating model skill at 
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reproducing extremal quantities requires some care. 
Clearly, the metrics used in traditional geostatistical 
analysis such as mean squared prediction error are 
unsatisfying for block-maximum data. Rather, we 
recommend the quantile score and the Brier score 
for threshold exceedences, as discussed and justified 
by Gneiting and Raftery (2007). These metrics are 
specifically tailored to evaluate the tail of the pre- 
dictive distribution and therefore seem more appro- 
priate in this context. 

2. TOWARD HIERARCHICAL BAYESIAN 
MAX-STABLE MODELS 

In their discussions of the relative merits of vari- 
ous approaches, the authors highlight the ability of 
hierarchical Bayesian models to represent richly flex- 
ible structures for underlying marginal parameters. 
As they point out, however, the conditional inde- 
pendence assumption made in the Bayesian analyses 
they discuss hamstrings the model's ability to pro- 
duce spatial association in process realizations. In- 
deed, others have also shown that failing to properly 
account for spatial dependence can lead to dramatic 
underestimation of uncertainty, and thus undercov- 
erage of posterior intervals, for the GEV parameters 
and return levels (Fuentes, Henry and Reich (2011)). 

The authors correspondingly laud the ability of 
max-stable processes to capture joint behavior across 
spatial locations, but lament the restriction to rel- 
atively simple underlying structures that pairwise 
likelihood fitting of max-stable process models im- 
poses. The trade-off between flexible marginal mod- 
eling and realistic spatial dependence modeling is 
almost treated as an inherent conundrum, almost 
analogous to a Heisenberg's uncertainty principle for 
spatial extremes. But we want to have it both ways! 

The authors rightly note that the unavailability of 
joint likelihoods for max-stable process models ap- 
pears to render their inclusion in hierarchical Bayes- 
ian models problematic. We view surmounting this 
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obstacle as a welcome challenge! As they note, prog- 
ress has already been made. For example, Ribatet, 
Cooley and Davison (2012) specifies such a hierar- 
chical model, but replaces the joint likelihood with 
a pairwise likelihood and modifies the resultant 
MCMC sampler using an asymptotic argument. The 
resultant sample from the "posterior" distribution 
appears to have desirable frequentist properties. 
While this approach may not be completely satis- 
fying in that it is computationally intensive and it 
does not pass the Bayesian purity test, it represents 
nice out-of-the-box thinking and is a clear step in 
the right direction. 

The authors also mention our recent manuscript, 
which describes an approach to hierarchical max- 
stable process modeling that we really like because it 
is fully Bayesian, straightforwardly produces predic- 
tions at unobserved locations, and can be fit to large 
data sets. Our approach suffers a bit because it does 
not readily generalize to most of the max-stable pro- 
cess models mentioned here by the authors. While 
possibilities certainly exist to expand on our ap- 
proach, and efforts are underway to do just that, 
we expect that completely novel angles and insights 
will be brought to bear on the rapidly-evolving field 
of Bayesian analysis for spatial extremes. 

3. SPATIAL MODELING OF HIGH 
QUANTILES 

Finally, we note that while the authors focus ex- 
clusively on the asymptotic extreme value thoery, 
other approaches do exist. Asymptotic arguments 
and parametric assumptions are clearly needed to 
estimate very extreme quantities, such as the 10,000 
year return level required by the Dutch Delta Com- 
mission. However, there are many important appli- 
cations where the focus is less extreme. For example, 
one may be interested in determining the effect of 
climate change on the 10-year return level of daily 
maximum temperature. Given the vast amount of 
meteorological data collected in the past century. 



there may be sufficient data to justify a less restric- 
tive and more robust approach such as quantile re- 
gression. Classical quantile regression (Koenker 
(2005)) gives a nonparametric estimate of covariate 
effects on a quantile (equivalent to a return level) of 
interest. Recently, semiparametric Bayesian quan- 
tile regression models have been proposed, includ- 
ing methods for spatial data (Lum (2010); Reich, 
Fuentes and Dunson (2011); Reich (2012); Tokdar 
and Kadane (2012)). An advantage of the Bayesian 
approach to quantile regression is the possibility of 
centering the prior of a flexible quantile model on 
a parametric extreme value distribution, and thus 
hopefully exploiting asymptotic arguments as the 
data deem appropriate. Comparing, and ideally merg- 
ing, quantile regression with extreme value analysis 
may be a promising line of future research. 
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